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Abstract. An extrinsic curvature surface model is investigated by Monte Carlo simulations 
on a disk. We found that the model undergoes a first-order transition separating the smooth 
phase from the collapsed phase. The results in this paper together with the previous ones suggest 
that the order of the transition is independent of whether the surface is compact (closed) or 
non-compact (open). 

1. Introduction This paper is aimed at showing that a surface model undergoes a 
first-order transition, which separates the smooth phase from the collapsed phase on a 
disk with an extrinsic curvature. The results in this paper together with previous ones 
[U [21 [31 m [S] indicate that the first-order transition occurs independent of whether the 
surface is closed or open. 

Two-dimensional surface emerges as an interface between two different materials in- 
cluding the air. Biological membranes are considered to be a two-dimensional surface sep- 
arating two-different biological materials [6l [3 [H [9l [10] , and they are described by some 
curvature Hamiltonian such as the one of Helfrich, Polyakov and Kleinert pTl [T2l [T3] . 

Fluctuations of surfaces are considered to soften/stiffen the surface [13], and theoret- 
ical investigations [Ill[ISl[I3[IHl[Tn] as well as numerical ones gD] [HI [H [H [H [25] [Ml 
\Tl\ [2ni [2H1 [301 [33 [33 [33] have been made on this problem. Surface models are classified 
into two classes: one is the extrinsic curvature model and the other is the intrinsic cur- 
vature model. It was reported that the intrinsic curvature models undergo a first-order 
transition [53 [13 [IS] , and we know that the transition is independent of whether the 
surface is closed [33 [28] or open [29] . 

The surface model with extrinsic curvature is also known to have a long-range order 
although the surface is softened by surface fluctuations. It was reported that the extrinsic 
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curvature model undergoes a first-order transition on triangulated spheres [TJ [U [31 3] as 
well as on triangulated tori 

However, it is still unknown whether the transition occurs on a disk with extrinsic 
curvatures. Many biological membranes are closed; the membranes play a role of an inter- 
face between two different materials, however, we know that there are holes in biological 
membranes. We must recall that a membrane with a hole is topologically equivalent to 
the disk. World surfaces of the string model can be divided not only into closed surfaces 
but also into open surfaces |10| . Moreover, it is still unclear how the open boundary 
influences the transition of the extrinsic curvature models. 

Therefore, it is interesting to study the phase structure of the surface model on a 
triangulated disk. In order to clarify this point, we study the surface model by using the 
Monte Carlo (MC) simulations; the model is defined by a Hamiltonian, which is a linear 
combination of the Gaussian bond potential and a bending energy. 




Figure 1. A triangulated disk of size (iV, Ne,Np) = (91, 240, 150), where iV, Ne,Np are the total 
number of vertices, the total number of edges, the total number of plaquettes (or triangles). 

2. Model Figure[l]is a triangulated disk of size A^ = 91, which is obtained by dividing 
each edge of the original hexagon into 5-pieces. By dividing each edge into i-pieces we 
have a lattice of {N,Ne,Np) = {3L^ + 31 + 1,91"^ + 3L,6L'^), where N,Ne,Np are the 
total number of vertices, the total number of edges, and the total number of plaquettes 
(or triangles). The flat and regularly triangulated disk such as the one shown in Fig[l]is 
used as a starting configuration of MC simulations. 

As stated at the end of the introduction we investigate the phase structure of an 
extrinsic curvature surface model on the triangulated disk. The model is defined by 
Hamiltonian S that is the linear combination of the Gaussian bond potential Si and the 
extrinsic curvature energy 

(1) S = XSi+bS2, Sl=Y,iX^-XJf , 

iv) 

i 3{i) 

where ^he sum over bond {ij) connecting the vertices Xi € R"^ and Xj E R'^, 

and A is the surface tension coefficient, which is assumed to A = 1, and b is the bending 
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rigidity. The vector n(i) in ^2 is the unit normal vector at the vertex i and is defined by 

(2) n(i) = N, = ^ ^m^i.) ' 

' Jii) 

where denotes the sum over triangles linked to the vertex i. The vector 

is the unit normal of the triangle and ^zij(i) is the area of the triangle We note 
that "^i includes the sum over the boundary bonds. As a consequence, the length 

of the boundary bonds can also vary in the simulations. Moreover, the boundary vertices 
share the bending energy 5*2 in Eq.(IT]), because 5*2 is defined not only on the internal 
vertices but also on the boundary vertices. 

In order to understand the interaction of 5*2 in Eq. ([1]), we graphically show in Fig. 
mja) the interaction range between the normals of triangles. Figure [2)Jb) is the interaction 
range defined by the conventional bending energy of the type 1 — • rij , where is the 
unit normal vector of the triangle i. The normal of the shaded triangle interacts with 
the normals of the other triangles in the figure. The range of interactions in Fig. Ufa) is 
relatively larger than that in Fig. [Hb). 




(a) The model in this paper (b) The conventional model 

Figure 2. The range of the interaction between the normal vectors in (a) the model of Eq.Q 
and in (b) the conventional model. The normal of the shaded triangle interacts with the normals 
of the other triangles shown in the figure. 

The partition function is defined by 

N 

(3) z{b) = / n ^^'^ ^'^p [-^(^)] ' ^(^) = ^1 + ^^2, 

where N is the total number of vertices as described above, and YiiLi ^-^i i^ ^ 3A^- 
dimensional integration. The expression S{X) shows that S explicitly depends on the 
variable X = {Xi,X2, • ■ • , X^), where Xi G is the position of the vertex i as described 
above. The center of the surface is fixed in the integration to remove the translational 
zero-mode. 

Note that the surface can be understood as the image of a mapping X from a two- 
dimensional parameter space P such as shown in FigUto R'^; X : P ^ R''. The surface 
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is allowed to self-intersect in our model and, hence it is called a phantom surface. This 
implies that the mapping X is not always injective. 

A physical observable (Q) is obtained by averaging the quantity Q{X) with the par- 
tition function such that 

N 

(4) ^^^=Z{b)J n dX^QiX) exp [~S{X)] . 

For Q — S2, we have the expression 

diogz _ idz 

where Z is dependent on the bending rigidity b and the surface tension A, which is 
assumed to be A = 1 as described above. The specific heat Cs^ can also be defined by a 
second- derivative of Z such that 

,r ^ b^d{S2) b^mogZ 

A phase transition is called first-order or discontinuous if the first-derivative of Z 
with respect to a parameter is a discontinuous function. Moreover, a phase transition is 
called second-order or continuous if the first-derivative of Z with respect to a parameter 
is continuous and the second-derivative of Z is discontinuous. Thus, the order of phase 
transition is originally connected to the analytic property of Z . 

We have another definition on the order of phase transitions. A phase transition can 
be called discontinuous if there is a physical quantity that discontinuously changes against 
a given parameter such as the temperature (or the bending rigidity in the surface model). 

Therefore, we can call the phase transition of the model as a first-order one not only 
if S2 is discontinuous against b but also if some physical quantity is discontinuous against 
b. 

3. Monte Carlo technique The 3iV-dimensional integration YiiLi '^-^i ^ '^^^ 
be performed numerically. The canonical Monte Carlo technique is used to update the 
variables X(e R^) so that X' — X + 5X, where the small change 6X is made at ran- 
dom in a small sphere in R'^. The new position X' is accepted with the probability 
Min[l, exp(— AS')], where AS — S'(new) — S'(old). The radius 6r of the small sphere is 
fixed at the beginning of the simulations to maintain about 50% acceptance rate. This 
procedure is called the Metropolis MC algorithm, and it produces equilibrium surface 
configurations X which satisfy the Boltzmann distributuion function eKjp[—S{X)]. 

We use a random number called Mersenne Twister [34] in the MC simulations. A 
sequence of random number is used for the 3-dimensional move of vertices X and for the 
Metropolis accept/reject of the update of X. 

By using a sequence of the equilibrium configurations of surfaces {X{i)} produced by 
the simulations, we have (Q) in Eq.Q simply by 

(7, (O^SfM^ 

In the following, we write (Q) as Q for simplicity. 
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We use the surfaces of size N = 25117, N = 15337, N = 9577, N = 5167, N = 2977, 
and N = 1387. The size of the surfaces is obviously larger than that of the same model 
on spherical surfaces used in [H [31 S] . 

The total number of Monte Carlo sweeps (MCS) after the thermalization MCS is 
about 4.7x10® close to the transition point of surfaces of size A^ = 25117 and iV=15337, 
3.0 X 10« - 2.6 X 10® for N = 9577 and N = 5167, and 1.5 x 10® ~ 1 x 10® for N = 2977 
and TV = 1387. Relatively smaller number of MCS is done at non-transition points in each 
surface. 

(a) 5 = 0.458 (collapsed) (b) 6 = 0.458 (smooth) 



(c) The surface section (d) The surface section 

Figure 3. Snapshots of the surface of A'' = 25117 obtained at 6 = 0.458; (a) a collapsed surface, 
(b) a smooth surface, (c) the surface section of (a), and (d) the surface section of (b). 

4. Simulation Results First, we show snapshots of surfaces of size N = 25117 in 
FigsEJa), and [3Jb) obtained at 6 = 0.458. The surface sections of them are shown in 
FigslSlJc), andlSUd). The surface in FiglHJa) is a collapsed surface, and that in FiglH^b) is 
a smooth one. 

FigurelDJa) shows the bending energy S2/2NB against b, where Nb is the total number 
of bonds including the boundary bonds. Dividing ^2 by 2Nb we have the bending energy 
per bond, where 2Nb comes from the definition of ^2 in Eq.([T]). We find that S2/'2Nb 
rapidly changes against b when N becomes large. This indicates that a phase transition 
occurs between the smooth phase and some non-smooth phase, although no discontinuous 
change can be seen in S2/2Nb- The reason why we can see no discontinuity in S2/2NB 
is because of the finitcness of N in the numerical simulations. 
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The mean square size is defined by 



(8) 



N ^ 



X,. - X) 



x^-Tx,, 

N ^ ' 



where X is the center of the surface, and are plotted in FiglDJb). No discontinuous 
change can be seen in in the figure, however, we find that the surface size rapidly 
changes as N increases. This reflects that the surface size discontinuously changes at the 
transition point in the limit of oo; the model undergoes the collapsing transition. 
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Figure 5. (a) The specific heat Cs2 vs. b, and (b) the peak values Csif^ vs. A'^ in a log-log scale. 
The solid line is drawn by fitting the largest three data to Ea. (|10|) and indicate a discontinuous 
transition. The dashed line is drawn by fitting the smallest four data to Ea. (|10p and is consistent 
with a continuous transition when TV < 9577. The error bars denote the statistical errors. 



The specific heat of the bending energy S2 can also be written as 

(9) Cs, =^((52-(^2))'), 
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and are shown in Fig. [HJa) against b. The errors shown as error bars in the figure are 
the statistical ones, which were obtained by the standard binning analysys. Anomalous 
behavior is obviously seen in at each N, and it indicates the existence of the surface 
fluctuation transition, since the anomalous peak of implies that the fluctuation of 
the surface becomes very large at the peak point. 

We plot the peak value Cg^^^ against N in FiglSlJb) in a log-log scale in order to clarify 
the order of the transition. The solid straight line in Figl5jb) was drawn by fitting the 
largest three Cg"^^ to 

(10) eg'''' ~ N", 
where cr is a critical exponent of the transition. We have 

(11) (T = 0.88 ±0.06, 

which is close to cr = 1. This implies that 5*2 discontinuously changes against b. Thus, 
the result indicates that the model undergoes a first-order surface fluctuation transition 
between the smooth phase and the non-smooth phase, where "the smooth" phase is 
characterized by small S2/2NB, and "the non-smooth" phase is characterized by large 
S2/'2Nb- The dashed line in FiglSlJb) is drawn by fitting the smallest four data to Eq. (fTO|) 
and is consistent with a continuous transition when N < 9577. This implies that the 
surface fluctuation transition is strongly influenced by the size effect. 
The collapsing transition can be characterized by the quantity 

(12) Cx^^^{{X^-{X')f), 

which is the variance (or the fluctuation) of X^. FiglBJa) shows Cx^ against b. Anomalous 
peaks can be seen in the figure, and they indicate a collapsing transition, which separates 
the smooth phase from the collapsed phase, where "the smooth" phase is characterized 
by large X'^ and "the collapsed" phase is characterized by small X^. 

Cx2 
0.2 



0.1 



0.455 0.46 0.465 iq^ ^ lo^ 

(a) b (b) 

Figure 6. The fluctuation Cx^ vs. b, obtained on surfaces of (a) — 25117, A'^ — 15337, and 
A'' = 9577, (b) a log-log plot of C^J" against A*' including the data obtained on the surfaces of 
Af = 5167, AT = 2977, and TV = 1387. 

Figure [6l^b) is a log-log plot of the peak value against N. The straight line was 
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drawn by fitting the data to 
(13) 

and then we have 
(14) 



V = 0.97 ±0.07. 



Thus, we confirni that the coUapsing transition is of first order. The exponent ~ 1 in 
Eq. (fT4| imphes that X"^ discontinuously changes against h. We note that the coUaps- 
ing transition is not influenced by the size effect in contrast to the surface fluctuation 
transition. 

S|/N 



1.502 
1.5 
1.498 



0:N=25117 
A:N=15337 
v:N=9577 



0.455 



0.46 , 0.465 
b 



Figure 7. The Gaussian bond potential Si/N vs. b obtained on surfaces of A'' = 25117 
Tfie expected relation Si/A'' = 1.5 is satisfied. 



9577. 



The Gaussian bond potential Si/N are plotted in Fig[7] in order to see that the 
expected relation Si/N = 1.5 is satisfied. The results clearly satisfy the relation, and 
therefore we consider that the canonical MC simulations were successfully performed. 



5. Summary and conclusion We have studied an extrinsic curvature model on a 
disk in order to see whether or not the model undergoes a first-order transition between 
the smooth phase and the collapsed phase. Monte Carlo simulations were performed on 
lattices of size iV = 25117 ~ 1387. 

We found that two types of transitions occur at the same transition point; one is the 
surface fiuctuation transition and the other is the collapsing transition. These transitions 
are well-known to occur in the same model on closed surfaces such as spheres and tori. 

The surface fiuctuation transition is strongly infiuenced by the size effect, and the 
simulation results obtained on the lattices of size less than TV = 9577 indicate that the 
model undergoes a continuous transition. However, the results obtained on the lattices 
larger than N = 9577 indicate that the transition is discontinuous. The peak height of 
the specific heat scales according to Cg^^^ ~ N°', and we have ct = 0.88 ± 0.06, which is a 
sign of the discontinuous transition of surface fluctuation. 

The collapsing transition is not influenced by the size effect in contrast to the surface 
fiuctuation transition. We have also obtained the fiuctuation Cx^ of X^, and find the 
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scaling property such that C™2^ . Thus we have j/ = 0.97 ± 0.07, which is a sign of 

the discontinuous transition of surface collapsing phenomena. 

Combining the results in this paper and the previous ones obtained on closed surfaces 
[U [51 [31 131 Ej J we conclude that first-order transitions can be seen in extrinsic curvature 
surfaces, and the order of the transitions is independent of whether the surface is closed 
or open. 



References 

[1] J-P. Kownacki and H. T. Diep, First-order Transition of Tethered Membranes in Three- 
dimensional Space, Phys. Rev. E 66 (2002) 066105(1 - 5). 

[2] H. Koibuchi, N. Kusano, A. Nidaira, K. Suzuki, and M. Yamada, First-order phase transi- 
tion of fixed connectivity surfaces, Phys. Rev. E 69, (2004) 066139(1 - 6). 

[3] H. Koibuchi and T. Kuwahata, First-order Phase Transition in the Tethered Surface Model 
on a Sphere, Phys. Rev. E 72 (2005) 026124(1 - 6). 

[4] I. Endo and H. Koibuchi, First-order Phase Transition of the Tethered Membrane Model on 
Spherical Surfaces, Nucl. Phys. B 732 [FS] (2006) 426 - 443. 

[5] H. Koibuchi, Phase transition of an extrinsic curvature model on tori, Phys. Lett. A 358, 
(2006) 339. 

[6] D. Nelson, Statistical Mechanics of Membranes and Surfaces, Second Edition, D. Nelson, 

T. Piran, and S. Weinberg Eds. (World Scientific, 2004) p.l. 
[7] G. Gompper and M. Schick, Self-assembling amphiphilic systems. In Phase Transitions and 

Critical Phenomena 16, C. Domb and J.L. Lebowitz, Eds. (Academic Press, 1994) p.l. 
[8] M.J. Bowick and A. Travesset, The Statistical Mechanics of Membranes, Phys. Rept. 344, 

(2001) 255. 

[9] K.J. Wiese, Phase Transitions and Critical Phenomena 19, C. Domb and J.L. Lebowitz 
Eds. (Academic Press, 2000) p. 253. 
[10] J.F. Wheater, Random surfaces: from polymer membranes to strings, J. Phys. A;Math.Gen 
27, (1994) 3323. 

[11] W. Helfrich, Elastic Properties of Lipid Bilayers: Theory and Possible Experiments, Z. 

Naturforsch 28c, (1973) 693. 
[12] A.M. Polyakov, Fine Structure of Strings, Nucl. Phys. B 268, (1986) 406. 
[13] H. Kleinert, The Membrane Properties of Condensing Strings, Phys. Lett. B 174, (1986) 

335. 

[14] D. Nelson, Statistical Mechanics of Membranes and Surfaces, Second Edition, D. Nelson, 

T. Piran, and S. Weinberg Eds. (World Scientific, 2004) p. 131. 
[15] F. David, and E. Guitter, Crumpling Transition in Elastic Membranes, Europhys. Lett 5(8), 

(1988) 709. 

[16] L. Peliti and S. Leibler, Effects of Thermal Fluctuations on Systems with Small Surface 

Tension, Phys. Rev. Lett. 54, (1985) 1690. 
[17] M.E.S. Borelli, H. Kleinert, and A.M.J Schakel, Quantum Statistical Mechanics of Nonrela- 

tivistic Membranes: Crumpling Transition at Finite Temperature, Phys. Lett. A 267, (2000) 

201. 



10 



T. ENDO ET. AL. 



[18] M.E.S. BorcUi and H. Kleinert, Phases of a stack of membranes in a large number of 
dimensions of configuration space, Phys. Rev. B 63, (2001) 205414. 

[19] H. Kleinert, Spiky Phases of Smooth Membranes. Implications for Smooth Strings , Euro. 
Phys. ,]. B 9, (1999) 651. 

[20] M. Bowick, S. Catterall, M. Falcioni, G. Thorleifsson, and K. Anagnostopoulos, The Flat 
Phase of Crystalline Membranes, J. Phys. I Prance 6, (1996) 1321; The Phase Diagram of 
Crystalline Surfaces, Nucl. Phys. Proc. Suppl. 47, (1996) 838; The Flat Phase of Fixed- 
Connectivity Membranes, Nucl. Phys. Proc. Suppl. 53, (1997) 746. 

[21] J.F. Wheater, The Critical Exponents of Crystalline Random Surfaces, Nucl. Phys. B 458, 
(1996) 671. 

[22] Y. Kantor, and D.R. Nelson, Phase Transitions in Flexible Polymeric Surfaces, Phys. Rev. 
A 36, (1987) 4020. 

[23] G. Gompper, and D.M. KroU, Phase diagram and scaling behavior of fluid vesicles, Phys. 
Rev. E 51, (1995) 514. 

[24] H. Koibuchi, A. Nidaira, T. Morita, and K. Suzuki, Monte Carlo Simulations of Branched 

Polymer Surfaces without Bending Elasticity, Phys. Rev. E 68, (2003) 011804. 
[25] H. Koibuchi, Z. Sasaki, and K. Shinohara, Phase transition of a tethered surface model with 

a deficit angle term, Phys. Rev. E 70, (2004) 066144. 
[26] H. Koibuchi, N. Kusano, A. Nidaira, K. Suzuki, and M. Yamada, Grand Canonical Monte 

Carlo Simulations of Elastic Membranes with Fluidity, Phys. Lett. A 319, (2003) 44. 
[27] H. Koibuchi, N. Kusano, A. Nidaira, Z. Sasaki and K. Suzuki, Phase transition of surface 

models with intrinsic curvature, Eur. Phys. J. B 42, (2004) 561. 
[28] I. Endo and H. Koibuchi, Phase transitions of a tethered membrane model on a torus with 

intrinsic curvature, Phys. Lett. A 350, (2006) 11. 
[29] M. Igawa, H. Koibuchi, and M. Yamada, Monte Carlo simulations of a tethered membrane 

model on a disk with intrinsic curvature, Phys. Lett. A 338, (2005) 433. 
[30] Y. Kantor, Statistical Mechanics of Membranes and Surfaces, Second Edition, D. Nelson, 

T. Piran, and S. Weinberg Eds. (World Scientific, 2004) p.lll. 
[31] Y. Kantor, M. Karder, and D.R. Nelson, Phys. Rev. A 35, (1987) 3056. 
[32] M. Bowick, A. Cacciuto, G. Thorleifsson, and A. Travesset, Eur. Phys. J. E 5, (2001) 149. 
[33] M. J. Bowick, Statistical Mechanics of Membranes and Surfaces, Second Edition, D. Nelson, 

T. Piran, and S. Weinberg Eds. (World Scientific, 2004) p. 323. 
[34] M. Matsumoto and T. Nishimura, "Mersenne Twister: A 623-dimensionally equidistributed 

uniform pseudorandom number generator" , ACM Trans, on Modeling and Computer Sim- 
ulation Vol. 8, No. 1, January pp.3-30 (1998). 



